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ABSTRACT 

The CI 1604 supercluster at z ~ 0.9 is known to contain at least four distinct member clusters, 
separated in both projection and redshift. In this paper we present deep, multicolor wide-field imaging 
of a region spanning ~ 45' on a side, corresponding to 21 h? Mpc (physical) at the supercluster 
redshift. We select galaxies whose colors correspond to those of spectroscopically confirmed cluster 
members in the i' versus r' — i 1 color-magnitude diagram. Using an adaptive kernel, we generate a 
map of the projected red galaxy density and identify numerous new candidate clusters that are likely 
supercluster members. Assuming that all of the density peaks are associated with the supercluster, 
its transverse size is ~ 10 h^Q Mpc, which is still significantly smaller than the nearly 93 h^o Mpc 
depth in redshift space. 

Subject headings: catalogues - surveys - galaxies: clusters: general - large-scale structure of the 
Universe 



1. INTRODUCTION 

The CI 1604 supercluster at z ~ 0.9 is one of the best 
studied high-redshift large scale structures. It was ini- 
tially det ected as two s epara te clusters in the plate-based 
survey of|G unn et all l)1986[) . with redshi fts and prelim- 
inary velocit y dispersions measured by iPostman et alJ 
(1998, 120011) . The proximity of the two clusters CI 
1604+4304 and CI 1604+4321 in both radial velocity 
(4300 km s" 1 ) and position on the sky (17') prompted 
iLubin et alJ ll2000j) to perfor m deep multi-band imaging 
with COSMIC l|Kells et al.lll99c| on the Palomar 5-m 
telescope covering the area between the two clusters. The 
imaging revealed additional overdensities of red galax- 
ies whose colors were consistent with spectroscopically- 
confirmed, early-typ e galaxies in the tw o original clusters 
(see Figs. 2 and 4 of ILubin et a,l.ll2?KKl) . suggesting that 
the clusters were part of a high-redshift supercluster. A 
total of four distinct red galaxy density pe aks were iden- 
tified in a contiguous area of 10'.4 x 18'.2. iGal fc Lubinl 
(2004) conducted a large spectrogra phic survey using 
LRIS and DEIMOS l|Faber et all2003fl on the Keck 10-m 
telescopes to confirm the new cluster candidates and pro- 
vide accurate velocity dispersions for all four supercluster 
components. They detected 230 total cluster members, 
confirming that all four components were in fact super- 
cluster members. These four clusters were nevertheless 
well-separated both in projection and in redshift, and 
have velocity dispersions ranging from 480 to 980 km 
s _1 , corresponding to Abell richness classes of R ~ — 2. 

The large radial depth of the supercluster (~ 93 
Mpc) in comparison to the small area imaged with COS- 
MIC (4.8 x 8.5 h^Q Mpc) suggested that this structure 
could extend to a significantly larger area on the sky. 
To test this hypothesis and search for new superclus- 
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ter components, we undertook a wide area, multi-band 
imaging survey usi ng the Large Format Camera (LFC; 
iSimcoe et al J 12000) on the Palomar 5-m telescope. In 
this paper, we present the results of this survey, includ- 
ing several new candidate supercluster members. We as- 
sume a ACDM cosmology with D, m — 0.3, A = 0.7, and 
H = 70 h 70 km s" 1 Mpc" 1 . 

2. THE IMAGING SURVEY 

Two fields in the region of the CI 1604 supercluster 
were imaged using the LFC on the Palomar 5-m tele- 
scope. The LFC is a mosaic of six 2048 x 4096 CCDs in a 
cross-shaped layout, mounted at prime focus. The imag- 
ing area corresponds roughly to the unvignetted region 
of a circle 24' in diameter. Data were taken on UT 2001 
May 15 and 17 and UT 2004 April 26 and 27, with the 
camera installed in the standard north-south orientation 
and all six CCDs in use. The imaging data were taken 
in unbinned mode, with a pixel scale of 0'.'182 pixel -1 . 
Seeing on all nights was better than 1'.'5, with average 
seeing of l'.'O. Individual exposures of 450s were taken, 
moving the telescope by up to 45" in each direction to 
alleviate the image gaps due to the spaces between the 
CCDs. Data were taken at two central pointings, using 
the Sloan Digital Sky Survey (SDSS) r',i', and z' filters. 
Only the night of UT 2004 April 27 was photometric, 
so individual exposures were taken in all three filters at 
both pointings to allow ca libration of the dee p images. 
A variety of standards from Smi th et al.1 l)2002fl were also 
observed over a range of airmasses on this night. Tabled 
provides the total integration times in each filter for the 
two pointings. 

Data reduct ion was perf ormed using the IRAF data re- 
duction suite (|Todvlll986|) . including the external pack- 
age mscred. 1 We followed the general mosaic imaging 

1 The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observatory, which is 
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TABLE 1 
LFC Imaging Log 



Pointing 



Coordinates, J2000.0 
RA Dec 



Total Integration (sec.) 
r' i' zf 



16:03:53.5 
16:05:05.7 



+43:21:33.3 
+43:04:18.4 



3600 
4950 



5850 
5850 



3600 
3600 



;uidelines used by t he NOAO Deep Wide Field Survey 
Jannuzi et alJl2004j) . with modifications appropriate to 
our data. 2 The data were overscan-subtracted, bias cor- 
rected, flat fielded, and fringe corrected in the standard 
way. Dark sky flats were generated by combining the 
deep images; the application of this secondary flat field 
results in images with less than 0.5% residual sensitivity 
variations across the field. Cosmic rays were detected 
using the craverage task, and a visual inspection of each 
image was performed to identify any additional bad pix- 
els or cosmic rays. All bad pixels due to cosmic rays, 
instrument defects, saturation, bleed trails, and satel- 
lite trails were masked. Astrometry was performed on 
each image using the task msccmatch to match detected 
objects with the USNO-A2 catalog. Typically between 
100 and 300 USNO-A2 stars are matched in each field 
and fitted with fourth-order polynomials in both dimen- 
sions, resulting in rms errors of ~ 0'.'3 in each coordinate. 
Individual exposures were then projected onto the tan- 
gent plane using mscimage, and all exposures of a single 
pointing in a single filter were combined with mscstack. 
Additionally, all exposures of each field (in r' , i', and z') 
were combined with mscstack to produce an extremely 
deep detection image. Images of photometric standards 
were subjected to the identical reduction as target expo- 
sures. 

Photometric calibrations were derived from standard 
star magnitudes measure d in a variety of circular aper- 
tures using SExtractor ijBertin k. Arnoutsl fl996) . We 
solved for zero points, as well as color and airmass terms 
using: 



+ A + Bsec(z) + C x {color) (1) 



where A is the zeropoint, B is the airmass coefficient, 
sec(z) is the airmass at which the standard was observed, 
C is the color coefficient, and the color used is r' — i' for 
the r' and %' data, and i' — z' for the z' data. 

Object detection in the target frames was performed 
by running SExtractor v2.3 in dual-image mode, using 
the ultra-deep image for detection, while measurements 
were performed on the single-band images. Photometry 
apertures are determined from the deep image, and these 
same apertures are used for the individual filter images. 
We use variable diameter elliptical apertures wit h ma- 
jor axis radius 2rfc, where Th is the Kron radius ijKronl 
fl980tlBertin" fc Arn outs 199(j); the magnitudes measured 
in these apertures are output as MAG_AUTO in SExtrac- 
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Fig. 1. — Positions of the two LFC and two COSMIC pointings 
in the CI 1604 region. 



tor. This procedure results in catalogs in which magni- 
tudes are measured using identical apertures in all three 
filters, which improves the measurement of galaxy colors 
by u sing the same phy sical size for each galaxy in all fil- 
ters lILubin et alJ20 00'l . To calibrate photometrically the 
deep exposures, shallower images taken on photometric 
nights were SExtractor-ed, and matched to the deeper 
data. Zero point offsets (typically ~ 0.1 mag) were de- 
rived and applied to the deep catalogs. To remove spu- 
rious objects around bright stars or near detector edges, 
as well as those with uncertain photometry, we require 
detection in both the r' and i' images, and exclude ob- 
jects with large photometric errors (MAG_ERR> 0.2 in 
either band), abnormally large semi-major axis (A > 10; 
typically faint objects deblended from a very bright com- 
panion), and zero isophotal area (objects in masked ar- 
eas). 

Because the two LFC pointings overlap slightly, the 
catalogs were combined and objects detected in both 
pointings had one detection removed. Additionally, there 
is a significant gap in the corner where the two LFC 
pointings meet. To fill this area, we used existing data 
from the C OSMIC imager, tak en as part of the survey 
presented in lLubin et alJ l)2000|) . Figure[I]shows the loca- 
tions of the LFC and COSMIC pointings. The COSMIC 
images were taken using the Cousins R and Gunn i fil- 
ters; to transform these data to the SDSS system, we 
matched objects in the overlap region between the COS- 
MIC and LFC fields, fitting relations for both the r' and 
i' filters of the form 



r' = A r R + B r (R-i) + C r 



(2) 



operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science 
Foundation. 

2 An extremely detailed data reduction description can be found 
online at the Palomar Observatory instrumentation website, at 
http: / /www. astro. caltech.edu/observatories/palomar/200inch/instruments. html 



The transformed COSMIC magnitudes show a la scat- 
ter of 0.1 mag relative to the LFC magnitudes. The cat- 
alogs from the two LFC pointings and the two COSMIC 
pointings were combined into a single master catalog, 
deleting duplicate detections in overlap regions with a 
preference for the slightly deeper LFC data. The final 
combined catalog contains ~ 25, 000 unique objects. 



3. CLUSTER DETECTION 
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Fig. 2. — The i' versus r' — i' color-magnitude diagram of all 
objects in our imaging area. Spectroscopically confirmed cluster 
members are overplotted as large filled circles. The vertical solid 
line shows the magnitude limit corresponding to the survey depth 
of r' = 24.5, i' = 24.5. The horizontal and vertical dashed lines 
show the color and magnitude cuts used to create the density 
maps in Figure The right panel shows the color distribution 
of all objects (dotted histogram) and spectroscopic members (solid 
histogram) scaled to the same peak number. 



The observed colors of early-type galaxies are strongly 
redshift dependent. Even at z ~ 1.0, the red sequence 
ijVisvanathan fc Sandaed 19771) of cluster ellipticals ca n 
be detected (iLuTfiTTeTaL 2000; IGladders & ^120001) . 
Thus, searching for overdensities of galaxies with the ap- 
propriate colors provides a simple yet effective means of 
detecting clusters even at moderate to high redshifts. In 
addition, we have obtained a su bstantial number o f red- 
shifts of supercluster members (|Qal fc Lubidl27)f7i) that 
can be used to verify the cluster galaxy colors. 

Figure plots the r' versus i' color-magnitude dia- 
gram (CMD) of all objects in our imaging area. Spec- 
troscopically confirmed cluster members are superposed 
as large filled circles. The vertical solid line shows the 
magnitude limit corresponding to the survey depth of 
r' = 24.5, i' = 24.5. The horizontal and vertical dashed 
lines show the color and magnitude cuts used to create 
the density maps in Figure 3. The right panel shows the 
color distribution of all objects with i' < 24.5 (dotted 
histogram) and the spectroscopically confirmed super- 
cluster members (solid histogram), scaled to the same 
peak counts. The red sequence of supercluster members 
is clearly visible at r' — i' ~ 1.1. Numerous bluer cluster 
galaxies are also evident; these are likely galaxies with 
ongoing star formation and therefore show strong [O II] 
A3727 emission redshifted into the r' filter. This hypoth- 
esis is consistent with the large fraction (~ 85%) of spec- 
troscopic members which show [O II] emission in their 
spectra. To perform cluster detection, we use only those 
objects with 1.0 < r' - i' < 1.4 and 20.5 < i' < 23.5. 
The color range was chosen to maximize the contrast 
of possible supercluster members relative to fore- and 
background galaxies, as seen in the right panel of Fig- 
ure |21 The i' magnitude limit was selected so that both 
the COSMIC and LFC imaging catalogs are reasonably 
complete in both r' and i' . These cuts eliminate the vast 
majority of foreground galaxies; only ~ 1500 objects sat- 
isfy these criteria. 

Using the color-selected galaxies, we produce an 
adaptive kerne l surface density map ({Silverman! Il986t 
iGal et al.l2003|) , with 10" pixels and an initial smoothing 
window of 128", corresponding to 1 hj^ Mpc at the su- 
percluster redshift. The adaptive kernel uses a two-stage 



process to generate a density map. It first produces an 
initial estimate of the galaxy density at each point in 
the map using a fixed grid. It then uses this pilot es- 
timate to apply a smoothing kernel whose size changes 
as a function of local density, with a smaller kernel at 
higher density. The smoothing size is chosen to sepa- 
rate the individual clusters in the supercluster without 
resolving much of the substructure in the clusters. The 
left panel of Figure[3]shows the density map for the color- 
selected objects from the combined LFC and COSMIC 
data, covering an area of 39'.8 x 43'.6. The four previously 
detected clusters at z ~ 0.9 are marked wi th circles and 
labele d with the same lettering scheme as IGal fc Lubinl 
( 2004). These four clusters are clearly detected as sig- 
nificant density peaks. However, there are a number of 
additional density enhancements of similar significance 
to the known members, especially to the northwest of 
the originally imaged region. To confirm that these struc- 
tures are not spurious, we generated several density maps 
using a variety of smoothing kernels corresponding to 
0.5 — 1.5 hjQ Mpc, as well as somewhat narrower and 
broader color cuts. Similar structures were visible in all 
such maps. We note the existence of numerous less sig- 
nificant density peaks and filaments in this map. These 
may be filamentary structure and possible very-low-mass 
systems associated with the supercluster. 

The right panel of Figure |21 shows the same density 
map, but only pixels with red galaxy surface densities 
of > 14,000 galaxies deg -2 are shown. This threshold is 
set to include all of the spectroscopically confirmed clus- 
ters. At least four additional structures, labeled "E," 
"F," "G", and "H", are above this threshold. We note 
that this choice of threshold is subjective; as mentioned 
before, many other overdensities are apparent. However, 
the spectroscopic observations indicate that the poor- 
est of the four confirmed clusters has a velocity disper- 
sion of only 489 km s" 1 , typical of low mass clusters 
ijBahcall fc Cenlll993|) . Therefore, the chosen density 
contrast will select galaxy concentrations consistent with 
even rather poor clusters. We also ran SExtractor on the 
density map with detection parameters fine-tuned to de- 
tect the four known clusters. This automated detection 
procedure picks up the same density peaks as those found 
visually. 

It is unclear from this map whether E and F are unique 
clusters, or substructures related to the known members 
C and D, respectively. Unfortunately, the newly detected 
structures are outside the region for which we have ob- 
tained spectroscopy, so no conclusions about their inde- 
pendence can be drawn at this time. Similarly, G and 
H are very close in projection and may be components 
of a merging system. They are, however, clearly sepa- 
rated from the other supercluster members. Their near- 
est neighbor is over 4 Mpc distant. Future DEIMOS 
spectroscopic observations are planned to determine su- 
percluster membership for the newly discovered struc- 
tures. 

4. CLUSTER PROPERTIES 

From the multicolor imaging data we are able to ex- 
amine a limited number of cluster properties for each 
candidate. The richness is estimated by counting the 
number of red galaxies (1.0 < (r' — £') < 1.4) in the mag- 
nitude range where both the COSMIC and LFC data 
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Fig. 3. — Adaptive kernel density maps of galaxies meeting our color criterion in the CI 1604 region. Previously known clusters are 
marked with circles. New candidates are marked with rectangles. The left panel shows the overall structure of the region. The right panel 
shows only areas with > 14000 red galaxies dcg~ 2 . Each map covers a 39(8 X 43(6 area. 



are complete (20.5 < i' < 23.5). If all of the new can- 
didate clusters are at the same redshift as the known 
ones, then this method samples the same (albeit bright) 
portion of the luminosity function and gives the relative 
richnesses of the clusters, subject to variations in the 
galaxy colors among the clusters. The computed rich- 
nesses, along with the cluster coordinates, are provided in 
Tabic |21 These coordinates are measured from the peaks 
in the red galaxy density map, and are slightly offset 
from t he spectroscopic centers reported in lGal fc Lubinl 
l|2004jh The photometric and spectroscopic centers of 
clusters A, C, and D agree to within 1": the photomet- 
ric center of cluster B is ~ 30" (0.25 h?Q Mpc) south of 
the spectroscopic center, likely because this cluster ap- 
pears significantly elongated to the south The r edshifts 
and velocity dispersions from lGal fc Lubinl l)2004[) for the 
previously detected supercluster members are also given. 
We note that cluster C has comparable richness to three 
other previously known clusters, despite having a much 
lower velocity dispersion. This cluster was only partially 
imaged in the earlier survey and thus may not have com- 
plete spectroscopic coverage. 



We also construct CMDs for each cluster, using a cir- 
cular region with radius 1 Mpc centered on each 
density peak. These are shown in Figure 0Jalong with 
the combined CMD for all eight clusters. To the right of 
each CMD, we show the color distribution in each field as 
the solid histogram; this is simply the result of collaps- 
ing each CMD along the magnitude axis. For compari- 
son, the dotted histogram shows the color distribution in 
a relatively blank region of the imaging data, scaled to 
an area with radius of 1 h^ Mpc. The enhancement of 
red cluster galaxies is clearly visible in these individual 
CMDs and color distributions. 



5. DISCUSSION 

If the red galaxy overdensities detected in this imag- 
ing survey are truly members of the CI 1604 supercluster 
at z ~ 0.9, we have identified by far the largest known 
structure at redshifts approaching unity. The largest pro- 
jected separation (AG) is approximately 9.6 /iy Mpc, 
still much smaller than the apparent 93 h^Q Mpc depth. 
It appears that we have detected the southeast limit of 
the supercluster; no significant overdensities are visible 
within ~ 10 h^Q Mpc directly to the east of A. We cannot 
rule out the structure continuing to the north/northeast 
or southwest. Most of the highest density peaks, includ- 
ing the four known clusters, are surrounded by smaller, 
lower density enhancements, suggesting a filamentary na- 
ture that is typical of large scale structure in simulations 
ijEvrard et alJl2002^ . The large range of densities in this 
supercluster makes it an ideal laboratory to study the 
dependence of galaxy properties and their evolution on 
their physical environment. 

We note that the apparent 8:1 axial ratio of this su- 
percluster is co mparable to or ev e n less than that found 
in simulations. iShandarm" et alJ l)2004|) provide an ex- 
cellent discussion of shape measurements for filaments 
and voids in numerical simulations with ACDM cosmolo- 
gies, demonstrating that superclusters typically have 
non-trivial, highly filamentary topologies with order-of- 
magnitude or greater differences in the lengths of their di- 
mensions (see their Figures 14 and 15). Furthermore, any 
complex filament will be easier to detect if it is viewed in 
projection along its longest axis. If the angular separa- 
tion between the sufficiently dense regions of a superclus- 
ter becomes larger than the contiguously imaged area in 
a survey, then it will only be detected as a single cluster. 

Our results demonstrate the efficiency of deep, moder- 
ately large area imaging in the vicinity of known clusters 
for detecting further structures. Hubble volume Simula- 
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TABLE 2 
Cluster Properties 



Cluster 


RA 


Dec 


Richness 


<Jl os (km s- 1 ) 


Redshift 


A 


16:04:20.9 


+43:04:40.6 


47 


962 


0.9001 


B 


16:04:23.6 


+43:14:30.5 


38 


719 


0.8652 


C 


16:04:06.5 


+43:15:40.4 


39 


489 


0.9350 


D 


16:04:35.7 


+43:21:19.6 


12 


640 


0.9212 


E 


16:04:05.0 


+43:13:31.9 


38 






F 


16:04:49.2 


+43:22:19.5 


37 






G 


16:03:45.4 


+43:24:07.0 


53 






H 


16:03:28.5 


+43:24:16.5 


58 










Fig. 4. — The i' versus r' — i' color-magnitude diagrams within a 1 h^ Mpc radius for the four previously known clusters and four new 
candidates, as well as the eight cluster areas combined. The dashed horizontal lines indicate the color range used for cluster detection. 
Histograms of r' — i' color for each cluster are plotted as the solid histogram to the right of each CMD, with the dotted histogram 
representing the color distribution in a blank sky region, scaled to the same area as the cluster histogram. 



tions show that similar superc lusters are not uncommon 
at z ~ 1 ijEvrard et al.11200^) . Therefore, we have un- 
dertaken a wide-field imaging survey of a large sample of 
z > 0.5 clusters to search for other large scale structures 
at these cosmologically significant redshifts. Using color 
cuts corresponding to early-type galaxies at a variety of 
redshifts we can also detect clusters at distances other 
than that of the target, i n a manner simila r to the cut- 
and-enhance technique of|Go^^^ayj2002 ) and the Red 
Sequence Cluster Survey l|Gladders fc YedboOOj) . 
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